FLACѧϰ2009-12-11 16:361. FLAC3DԪ
𣺲ǣ޲ַ

2. ҪյЩ
Ҫgen, ini, app, plo, solveȽģʼ߽

3. ģ͵ӣ
plo blo groԿͬgroupɫֲ

4. ģ͵ı߽
plo gpfix red

5. ģ͵ֲ
plo fap red

6. ģ͵ͼ
λƣplo con dis (xdis, ydis, zdis)
    Ӧplo con sz (sy, sx, sxy, syz, sxz)

7. ģ͵ʸͼ
plo dis (xdis, ydis, zdis)

8. ģжٵԪڵ㣿
pri info

9. ģ͵ĺͼ
File/Print type/Jpg fileȻѡFile/PrintʽѡΪjpeļ

10. һļ
File/callcall

10. ʩ
app nstress

11. εͼĴСǶȣ
ۺʹx, y, z, m, ShiftʹCtrl+RCtrl+Zȿݼ

12. νб߽Լ
fix x ran ԼٶȣڳʼԼЧλԼ

13. ֪ÿԪID
˫Ԫı棬֪ԪID

14. νƬ
plo set plane ori () norm (ʸ)
    plo con sz plane (ʾzӦƬ)
    
15. α
save +ļ.

16. εѱĽ
rest +ļFile / Restore

17. ͣ㣿
Esc

18. ڳнͣɻָ㣿
мpausecontinueм

19. ĳ㲽
ڼаոμ㣬Զһ

20. Fishʲô
FLAC3Dԣвģ͡ܽеĹ

21. FishǷһҪѧ
𣺿ԲãҪʱMannualҪıͿ

22. FLAC3DļʽЩ
νֻҪıļʲô׺

23. εһЩѡģ飿
config dyn (fluid, creep, cppudm)

24 .ԲʩԼ
𣺿fix ... ran cylinder end1 end2 radius r1 cylinder end1 end2 radius r2 not

r2<r1ʵѡrangeʱѡԲĲõ߽硣

25.ܰһPLOTͼݵԱͼ
set log on ݵתexcelﴦһ£Ȼsurferʲôͼ

ƾˡ

26.dddip⣺
 Ҫȷ㶨ķ߷dip ǵ˼xoyͶļнǣdipǣdip ʱ۽ǡddy˳ʱתxoyϵͶӰʱתĽǶȡ

27.εFISHļ
ֵFISHļԭǣļΪ fishcall.fis.txt(Ϊڼ±txtĺ

׺Ĭϵģ޸ķǣ fishcall.fis.txtȫƵµļ±ٰ

ʱ򽫱͸Ϊļ(*,*),ʱļ fishcall.fisҲ

ȷˡ

28.Tension-p shear-p tension-n shear-nʾĺ壿
ڼѭ,ÿѭ,ÿzonefailure criterionڲͬ״̬,shear 

tensionֱʾܼ״̬.nʾnow,pʾprevious,ֱʾںǰ

ѭд״̬.

29.һЩ
plot hist m vs nʽ,my,nx(m,n);
"-"ʾֵ"mirror",ʹAUTOCADĻ"mirror"ĺ

30.β鿴һλֵ
print gp disp range x() y( ) z( )
print gp disp range id 
plot gp disp range x() y( ) z( )
plot gp disp range id

31. a: fix x y z range x -0.1 0.1 
    b: fix x range x -0.1 0.1    aʽbʽк
a ʾ x=0 ƽ ǹ̶ģb ʾx=0 ƽ ֻxܷ˶

32.initial  apply к
𣺳ʼʼӦ״̬ȣ
apply߽ʩӱ߽λƵԼȡ
initialӦ״̵̬ķı䣬һҪʼapplyʩӵı߽

ᷢ仯

33.β鿴ʱβƽľֵ
𣺲his¼㣬λӦ his unbal 
his gp(zone) zdis range (0 0 0) id=
his write n vs m begin ʱend ʱ +file filename.his
Ҫƽľֵ
his unbal
step 100000 or solve
his write 1 vs step begin 1 end 1000 file 123.his
ʹͿԲ鿴ʱµĲƽľֵ

34.ʾߵ 
plo ske magf 10   10ΪŴϵ

35ʼţ
𣺳ѧϰȿ˵ User's guridSection2 Getting startedSection 3 Problem

sloving with FLAC

36ЩҪգ
grid gen prop interface ini fix set hist step solve pause cont sav ȣѧϰĹ

ЩͲӦò鿴˵Command Reference潲ĺرҪע

һҪϸһCommand ReferenceͷOptional Range PhrasesǽΧģΪ

ܶʩӶҪӵģ͵һ֣͸ʾһ֡һԭͬʺ

5ITASCAһFLAC FALC3D UDEC 3DEC PFC PFC3Drange

37FLACFLAC3Dģʲô
FLACFLAC3DĽģ˼һģͬFLACGRID񻭳

ȻMODELL NULLֳಿ֣GENЩһһȻ

GEN LINEGEN ARCҪڵı߽磬FLACвֻһGRIDֻ

һΡFLAC3DֱGEN һһӣGENúܶΡ

38鿴ģͣ
plot grid Բ鿴plot grid num Բ鿴ڵš

39鿴߽
plot fix

40鿴ԤĿڱ߽磿
plot grid mark

41ô鿴ģ͵ͼ
pl syy (sxx xdis ydis sig1..) fill ע⣬fillؼ־ǲ鿴ͼfill ؼ

ǲ鿴ֱͼ

42鿴ģ͵ʸͼ
pl grid vel (disp) red

43鿴INTERFACE
pl grid iface id red 
print iface id ; id interfaceı

44FLACôͼ
ctrl+z shift+z ,ȻҪŴλϽǺ½ǾͿԷŴĳλҪС

˫Ļ

45жģһķ
   ʾģIDIDţ˵ϵӦIDš

46ıģӦ߽ʱFLAC3DECĲ֮ͬ
3DECУıӦ߽ʱģͱӦֱӱFLACиıӦ߽ʱģ
ӦΪӦӡ

47flacӰʱ
1Ԫϣṹ֮ĸնȶԱȺܴ
2ֵԪߴԱȺܴ󣬾ⳤϸȴ51ϸԪ

48 FLACжвģ
11бģ%_4[ W/N,n R5q
1ģ    ڿڻ
2ͬԵģ
3Եģ
4۸ͬԵģ
5³ˣģ
6Īģ    һѧΪ
7ģ    оƻΪȶú󿪲ɣ
8ӦӲ/ģͣģ   ״ƻѧΪ
9˫ӦӲ/ģ ٽ״㿪
10˫ģ
11Ľճģ

49Բƽ
flacƽ͵ıֵRʾģ͵ĲƽԴСֵΪٷҴ

С㣬ˣģͲﵽȶƽ״̬ԸҪľ趨RֵģĬ

Rֵ1e-5
עRֵСʾƽ⣬ģͻܷͨ plot vel 鿴
ٶȣһ漴ٶȳûƽ״̬Ŀɿָʾ


50ͣ
plot add cont disp out on ʱʱʾͼ
plot block group gpnum on ʾڵ
plot inter red ske         ӴЧ.
del range group 1 not ɾgroupϢ    
del range group 2 ɾgroupϢ
ini z add ֵС range group name ˳zƶ

Ӵ interface 1 face range x(Χy(Χ z(Χ
           interface 1 face range plane orig(x y z ꣩ norm߷
water table clear ˮλ
pl set back white ;ñ
set plot bitmap size (800,510) ;ͼƬС 
plot set caption size 37 ; ͼִС (default=35 Χ10-50)
plot set center 100 0 40 
fix x y z range x -0.1 0.1 ʾx=0 ƽǹ̶ fixed boundary
fix x range x -0.1 0.1ʾx=0 ƽֻxܷ˶roller boundary

51fishĽͣ
def aux1 庯
loop n(1,nptab1) ʼһѭ1ѭnptab1
rr=xtable(1,n) table 1 xжnrr
zz=ytable(1,n) table 1 yжnzz
loop k(1,nprof) ʼڶѭ1ѭnprof
case_of k һѡṹñkж
i F alfa=0.5*(pi+.1)*float(k-2)/float(nprof-2)
xtable(n+10,k)=rr*cos(alfa)
ytable(n+10,k)=rr*sin(alfa)
case 1 k1ʱִ
xtable(n+10,k)=rr table n+10xеkֵڱrr
ytable(n+10,k)=-40.0 table n+10yеkֵڱ-40
case 2 k2ʱִ
xtable(n+10,k)=rr table n+10xеkֵڱrr
ytable(n+10,k)=0.0 table n+10yеkֵڱ0
end_case ѡṹ
end_loop   ڶѭ
end_loop һѭ
end 
aux1 ִʲôĺ


52. ڹ̶ʩλʱfixapplyͬʱʹãӦinifix,ʩϡ

53.˵ plot items/add/structure elements/geometry/ṹԪʽѶԻmag facֵĳɲͬͿԿλƷŴĽṹԪ

54.üϵϡ¡Ҽڴƶģ͡

55.ܰһPLOTͼݵԱͼ
set log on ݵתexcelﴦһ£Ȼsurferʲôͼƾˡ


56 ڽѹʱʩΧѹ
ʹapply sxx =-1e5 szz=-1e5 range cyl end1 0,0,0 end2 0,2,0 radius 1XZ۾Ե档

57 ڽѹʱζλƿƣ
ʹ ini yvel 1e-7 range y -.1 .1 ini yvel 1e-7 range y 1.9 2.1 ɣֵԱ仯

58 ڽб¼ʱΪʲôsol, fos вУ
Ӧsol fos slope3dfos.sav associated(ҪҲУ restore üɡ

59 FLAC3DУμ¼еֵʱʱ¼ôŪ 
 hist wtite n1 vs n2 file.txt ¼set hist_rep n

60 ¼Ӵϵѧô죿
hist inter i nstress /ndis/sstress/sdisp

61 Ӵòļ
coh dial fric kn ks ten.

62 Ԫݣ

expgrid .flacsd

63 ĳͼϵʸֵ D Q Z9s1} e4Y U%q
   plo add ske 
plo add dis plane
plo add axe

64 Ƭ
    plot set plane ori x y z norm x y z

65 һ֧ڵӣ
   gen zone radbrick p0 (0,0,0) p1 (10,0,0) p2 (0,10,0) p3 (0,0,10) &
size 3 5 5 7 rat 1 1 1 1.5 dim 1 4 2 fill
model mohr
prop bulk 1e8 shear .3e8 fric 35
prop coh 1e3 tens 1e3
fix x range x -0.1 0.1
fix z range z -0.1 0.1
fix y range y 9.9 10.1
set grav 0 10 0
ini dens 1000 
ini syy=0.0 grad 0,-10000,0
ini sxx=0.0 grad 0,-5000,0
ini szz=0.0 grad 0,-5000,0
apply sxx=0.0 grad 0,-5000,0 range x 9.9,10.1
apply szz=0.0 grad 0,-5000,0 range z 9.9,10.1
set large
hist unbal
solve
save tr_eq.sav
; monitor trench excavation
ini xdis=0 ydis=0 zdis=0
hist gp xdisp 1,0,0
hist gp zdisp 0,0,2 ;
; excavation steps 
range name exc1 x 0 1 y 0.0 0.8 z 0 2 
range name exc2 x 0 1 y 0.8 1.6 z 0 2
range name exc3 x 0 1 y 1.6 2.4 z 0 2
range name exc4 x 0 1 y 2.4 3.2 z 0 2
range name exc5 x 0 1 y 3.2 4.0 z 0 2
; excavation step 1
model null range exc1
step 10000
save exc1.sav
; excavation step 2
model null range exc2
step 1000
save exc2.sav
; excavation step 3
model null range exc3
step 1000
save exc3.sav 
; excavation step 4
model null range exc4
step 1000
save exc4.sav
; excavation step 5
model null range exc5
step 1000 
save exc5.sav
; install cable support
rest exc2.sav
sel cable beg 1.0,0.4,1.5 end 5.0,0.4,1.5 nseg 4
sel cable beg 1.0,0.4,0.5 end 5.0,0.4,0.5 nseg 4
sel cable beg 1.0,1.2,1.5 end 5.0,1.2,1.5 nseg 4
sel cable beg 1.0,1.2,0.5 end 5.0,1.2,0.5 nseg 4
sel cable prop emod 2e9 ytension 1e8 xcarea 1.0 gr_coh 1e10 gr_k 2e9 gr_per 1.0
model null range exc3
step 1000
save cab3.sav
sel cable beg 1.0,2.0,1.5 end 5.0,2.0,1.5 nseg 4
sel cable beg 1.0,2.0,0.5 end 5.0,2.0,0.5 nseg 4
sel cable prop emod 2e9 ytension 1e8 xcarea 1.0 gr_coh 1e10 gr_k 2e9 gr_per 1.0
model null range exc4;
step 1000
save cab4.sav
sel cable beg 1.0,2.8,1.5 end 5.0,2.8,1.5 nseg 4
sel cable beg 1.0,2.8,0.5 end 5.0,2.8,0.5 nseg 4
sel cable prop emod 2e9 ytension 1e8 xcarea 1.0 gr_coh 1e10 gr_k 2e9 gr_per 1.0
model null range exc5
step 1000
save cab5.sav
plot create TrenchSupport 
plot set rot 63.89, 334.35, 3.4
plot set cent 2.893, 3.446, 3.899
plot set size 4.136
plot add cont disp out on
plot add sel cable force
plot add axes black


66 һѹӣ
gen zone cyl p0 0 0 0 p1 1 0 0 p2 0 2 0 p3 0 0 1 size 4 5 4 
gen zone reflect norm 1,0,0
gen zone reflect norm 0,0,19
model mohr 
prop bulk 1.19e10 shear 1.1e10
prop coh 2.72e5 fric 44 ten 2e5
fix x y z range y -.1 .1
fix x y z range y 1.9 2.1
apply sxx =-1e5 szz=-1e5 range cyl end1 0,0,0 end2 0,2,0 radius 1
ini yvel 0 grad 0 -1e-4 0 range y 0 2
sol



1.FLACʽʹǰ׼

1

ֽϻصλáزϡ򵥱ʾ뼰ϡ

2ϳͬһλ

еزϣ Density, Bulk modulus, Young`s modulus, tension, cohesion, friction Angle ϣͬһ λ

ע ϵĵλλͬɹСֵ FLAC޷㣬ERRORϢ  

3ӦùʽԼ

Ӧùʽݾ飬FLACϵķΧΪDebugϷʱ֤

4x, ynode i, j֮Ĺϵ

ͼϾ뼰ȵĹϵx, yϵͳx, yϵͳתϵΪ˱ԺϵķӦȷʵ λüӦx, y ꡣ

ʹýܵطʹý񣬸տʼִгʱʹ̫Ŀʹ 1000Խʡʱ䡣

2.FLAC 3D ıд裺

1 Config ________

2 Grid ________

3 Model ________

4 ʼӦƽ________

    (1)x, yĹϵʹGenָʾ

Gen x1,y1 x2,y2 ,x3,y3 x4,y4 i=i0,i1 j=j0,j1 ϸָμʹֲᣬFLACʽԶx, y Ĺϵ ڲ겻׿ƣĹӰ죬ʹǰ4Ĺϵ̶ܳҪĲͬʵ ʿꡣ

    (2)趨ʣprop

    (3)趨Set Grav, Apply Pressure, ini sxx, Syy

    (4)趨߽fix, free

    (5)ʼӦƽ⣺solve

    (6)棺Save

5 󹤳̵Ӱ________

ڵӦ̵ֲ̼׶ִи̽й̵Ӱ죬ʹµĲ裺

    (1)Ӧƽ⣺re_____ .sav

    (2)趨µĲʣmodel,prop

    (3)趨µ֧ʣstruct

    (4)趨µ

    (5)趨߽

    (6)󹤳ʱӦƽ

    (7)
 
